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Abstract 


How to distinguish frauds from consequences of strategic behavior is the primary challenge 
for election forensics. Election forensics methods use empirical distributions of turnout and 
vote choices to identify fraudulent activity. Strategic behavior can affect these 
distributions. Generally, many types of election forensics methods interpret a multiplicity 
of modes in election data as indicators of frauds, but strategic behavior induces correlations 
among electors’ behavior that can produce multimodalities. We present agent-based models 
(ABMS) designed to represent electors who decide how to act in elections based both on 
their own tastes and on their beliefs about what the other electors will do. In particular we 
define ABMs designed to match equilibria derived in models of wasted vote logic and of 
strategic abstention. To facilitate computing pivotal vote probabilities, our ABMs use a 
Poisson games formulation of the models instead of the fixed-electorate assumptions the 


formal models originally used. We replicate the formal models’ basic equilibria. 


1 Introduction 


A key challenge for election forensics—the field devoted to using statistical methods to 
determine whether the results of an election accurately reflect the intentions of the 
electors—is to be able to distinguish election results caused by election frauds from results 
produced by strategic behavior or other normal politics (Mebane 20134, 2016). Election 
forensics studies counts of votes, counts of eligible voters and other traces of an 
election—preferably at low levels of aggregation such as tallies for each polling station—to 
produce evidence regarding what happened in the election. By starting with the numerical 
results and other measures of the voting process election forensics does not address the 
entirety of an election or address the full range of frauds that are possible (Lehoucq and 
Jiménez 2002; Lehoucq 2003; Magaloni 2006; Schedler 2006; Levitsky and Way 2010; 
Minnite 2010; Birch 2011; Hyde and Marinov 2012; Svolik 2012; Wang 2012; Simpser 2013; 
Stokes, Dunning, Nazareno and Brusco 2013; Norris 2014). For example, if parties are 
excluded from the ballot, such an action may not produce distinctive patterns in the votes 
that are cast. But some violations of electoral integrity such as unfair access to campaign 
resources, wrongfully manipulated voter lists, vote buying, voter intimidation and other 
coercive actions may produce distinctive patterns in votes that statistics can detect. 
Ambiguities arise when such patterns might also be produced by strategic voting and other 
normal political activites such as forming coalitions. 

Many methods for trying to detect election frauds have been proposed (e.g. Myagkov, 
Ordeshook and Shaikin 2009; Levin, Cohn, Ordeshook and Alvarez 2009; Shikano and 
Mack 2009; Mebane 2010; Breunig and Goerres 2011; Pericchi and Torres 2011; Cantu and 
Saiegh 2011; Deckert, Myagkov and Ordeshook 2011; Beber and Scacco 2012; Hicken and 
Mebane 2015; Montgomery, Olivella, Potter and Crisp 2015; Mebane 2016; Rozenas 2017; 
Ferrari and Mebane 2017; Ferrari, McAlister and Mebane 2018). Methods based on the 
second significant digits of vote counts have been shown to respond both to normal 


political activities (strategic behavior, district imbalances, special mobilizations, coalitions) 


and to frauds (Mebane 2013a, 2014). Methods that examine the last digit of vote counts 
can be fooled if malefactors have sufficient control over the numbers (Mebane 2013b). All 
the methods can effectively identify various kinds of anomalies, but assessing whether the 
anomalies are due to frauds presents further challenges. 

Some of the methods explicitly focus on the modality of election data. Some methods 
in this vein emphasize that unproblematic elections feature unimodal distributions of 
turnout and regular flows of votes—the latter are most compatible with assumed unimodal 
distributions for parties’ shares of the votes (Myagkov, Ordeshook and Shaikin 2008, 2009; 
Levin et al. 2009). Other contributions connect “spiky” (hence multimodal) distributions 
of turnout and vote proportions to ideas about agents committing frauds in ways that they 
intend to be detected (Kalinin and Mebane 2011; Mebane 20130; Rundlett and Svolik 
2015). The sharpest contribution featuring multimodality is a model proposed by Klimek, 
Yegorov, Hanel and Thurner (2012) that stipulates a particular functional form according 
to which frauds occur. 

Review of the theory (Borghesi 2009; Borghesi and Bouchaud 2010) that motivates the 
Klimek et al. (2012) conception suggests, however, that multimodal distributions may be as 
readily produced by strategic voting, coalitions and other strategic behavior as by election 
frauds that stem from maleficent activity (for details see Mebane (2016)). That theory, as 
it has been developed so far, does not imply that the distributions produced by strategic 
voting and by frauds are the same: currently the theory is not specifically quantitative. But 
the theoretical ambiguity about the origins of multimodality may carry over to make the 
parameters of the Klimek et al. (2012) conception ambiguous. Many of the other statistical 
methods also, on inspection, are crucially triggered by multimodalities and so may be 
subject to similar ambiguities. An election may appear to have frauds when in fact it has 
only robust politics featuring strategic activity by electors—by voters and would-be voters. 

Here we report on the initial steps of our plan to use agent-based models (ABMs) to 


produce quantitatively specific realizations of theories of strategic election behavior and of 


election frauds. We hope to use ABMs to simulate various election systems at realistic 
scales, with agents that are involved in various kinds of networks and endowed with various 
distributions of attributes, where agents engage in equilibrium behavior such as is studied 
in formal models of rational election behavior (e.g. Palfrey and Rosenthal 1985; Cox 1994). 
We will introduce various kinds of frauds into such systems. Among the questions will be 
to what extent can statistical methods for election forensics discriminate strategic actions 
from frauds, and can the methods measure the extent and location of the frauds accurately. 

Previous efforts to use simulations motivated by strategic voting ideas to develop data 
to use to assess election forensics methods have not rigorously imposed equilibrium 
conditions, but have only roughly tried to approximate expected consequences of such 
behavior. For instance, simulations have shifted votes between candidates with references 
to the kind of behavior that happens when voters act according to wasted vote logic (e.g. 
Mebane 2013a). Perhaps such approaches capture some important qualitative aspects of 
what happens when voters act that way, but there is no proof of that and there is little 
reason to believe the quantitative details about how many votes are affected are correct. 
While some formal models have equilibrium ideas with notions of frauds (e.g. Simpser 
2013; Rundlett and Svolik 2015), no simulations have attempted to achieve quantitative 
precision about what happens. 

Here we take only the first steps toward such a plan. We use ABMs to replicate 
equilibria obtained in two important formal models of elector behavior: the Palfrey and 
Rosenthal (1985) model of strategic abstention, as updated by Demichelis and Dhillon 
(2010) to include aspects of learning; and the Cox (1994) model of wasted vote behavior. 
In both instances we modify the models to use the conception of Poisson games (Myerson 
1998, 2000). Others have studied election equilibria using the Poisson games formulation 
(e.g. Bouton and Gratton 2015). We take a Poisson game approach to make it feasible to 


compute pivotal voting probabilities (Skellam 1946). 


Game Theory and Agent Based Models 


Social science is not just concerned with individual behavior, but also, with how individual 
interactions aggregate to larger scale events (Axelrod and Tesfatsion 2006). To this end, 
agent based models allow us to study systems with the following characteristics: “(1) the 
system is composed of interacting agents; and (2) the system exhibits emergent properties, 
that is, properties arising from the interactions of the agents that cannot be deduced 
simply by aggregating the properties of the agents” (Axelrod and Tesfatsion 2006). Indeed, 
the strength of ABMs is their ability to capture interdependent human behavior (Bruch 
and Atwell 2015; Epstein 1999). 

That is, individuals behave in a social context that is defined by the collective behavior 
of individuals. There is feedback between micro level behavior on the one hand, and macro 
level context on the other: “In the short run, individuals respond to their environments; in 
the longer run, the accumulation of individuals choices or behavior changes the 
environment” (Bruch and Atwell 2015). Both game theory and agent based models 
attempt to capture this interdependence; player strategy is contingent on other players’ 
actions and agent behavior is determined, in part, by interactions with other agents. 

It is worth noting the similarities and differences between game-theoretic models and 
ABMs. Both techniques assume an explicit payoff function, focus on micro-level behavior, 
and isolate specific components or parameters while holding details fixed (Johnson 1999; 
de Marchi and Page 2014). Computational models focus on dynamics, while game theory 
focuses on equlilibria (de Marchi and Page 2014). Both types of models constrain behavior, 
but while game theorists do so to make equilibria analytically tractable, computational 
models draw constraints from cognitive, or decision-making research (de Marchi and Page 
2014; Kim, Taber and Lodge 2010). Rational choice is included among the cognitive 
heuristics found in ABMs and as such, game theoretic and agent based models are not 
inherently opposed (Axelrod 1997; Johnson 1999; de Marchi and Page 2014). The 


immediate advantage of ABMs is the ability to move beyond the rational choice 


assumption to model behaviors that are analytically intractable (Axelrod 1997). 

Both techniques have mechanisms to study changes in parameter space. For game 
theory, this is comparative statics. For ABMs, we relax assumptions about micro level 
behavior to test if macro level outcomes occur as expected (Epstein 1999). This is 
invaluable for theory testing. Robustness in a computational model is important, as we 
must avoid “brittleness.” That is, equilibria should be invariant to changes to initial 
parameters (Holland 1983). Additionally, in replicating formal theoretic models in an ABM 
framework, we find that computational power complements mathematical precision, 
leveraging the comparative advantage of each technique to gain greater insight than either 
one alone can provide (Miller and Page 2004). 

This synthesis of technique—analytic on the one hand and agent-based on the 
other—parallels behavioral game theory, an approach that also has tried to solve the 
“paradox of voting” (Downs 1957; Bendor and Diermeier 2003). Previous research on 
elections employs adaptive rationality, whereby agents learn not by optimizing utility, but 
rather, by trial and error. This is implemented via reinforcement learning, where agents 
pursue actions that have previously given satisfaction and inhibit those that have proven 
unsatisfactory. (Kollman, Miller and Page 1992; Bush and Mosteller 1955; Bendor and 
Diermeier 2003). Note that previous work has focused on parties as agents—acting either 
ideological or office-seeking—to explore how polling affects party strategy (Kollman, Miller 
and Page 1992). This project focuses on voters as agents to explore how perceptions of 


electoral closeness affect turnout.! 


1.1 Strategic Abstention 


Palfrey and Rosenthal (1985) consider the classic instrumental formulation for the net 


‘We recognizes that some distinguish equation based models (EBMs) from ABMs (Parunak, Savit and 
Riolo 1998). EBMs begin with equations that capture relationships between observable characteristics (party 
preference, for instance), while ABMs begin with behaviors that define how individuals interact with one 
another (Parunak, Savit and Riolo 1998). The literature is divided on whether this differentiation makes 
sense; for any given computational model, there is a set of equations that could (theoretically) produce the 
same result (Bonabeau 2002; Epstein 1999). For the purposes of this paper, we do not differentiate. 


benefits R from voting in an election with two alternatives: 


R=pB-C+D, (1) 


where p is the probability that an individual’s vote decides the election outcome (is 
“pivotal” ), B is the difference in benefits for the individual if the individual’s more 
preferred choice wins instead of the alternative, Cis the cost of voting and D is the direct 
benefit from voting. Neither C nor D depend on the election outcome but depend only on 
the individual deciding to participate in the election by voting. Palfrey and Rosenthal 
(1985) demonstrate the existence of equilibria with positive turnout given a variety of 
assumptions about the distributions of costs and uncertainty. Palfrey and Rosenthal (1985) 
argue that turnout, when positive, is usually low. 

Demichelis and Dhillon (2010) consider symmetric equilibria given the same net 
benefits and structure as Palfrey and Rosenthal (1985), while adding learning dynamics. 
Demichelis and Dhillon (2010) normalize (1) by setting B = 1, in which case a necesssary 
condition for an individual to vote is p — 2c > 0, where c = C — D. Letting gq denote the 
symmetric mixed strategy for all players and using g(q, M) to denote the probability of 


being pivotal, Demichelis and Dhillon (2010) suppose that qg changes over time according to 


apo hw) (2) 


where sign K(q) = sign(g(q, M) — 2c). Demichelis and Dhillon (2010) investigate the 


dynamic stability of equilibria. 


1.2 Wasted Vote Logic 


Cox (1994) considers situations where each voter casts a single vote, M alternatives “win” 
the election by gaining seats and the M alternatives that have the M highest numbers of 


votes win the seats. Assuming that each voter maximizes expected utility based on pivot 


probabilities, that utilities are diverse,” that the distribution of preferences is common 
knowledge and that expectations about the election outcome are publicly generated (to the 
point that Cox (1994, 610) imposes a rational expectations condition), Cox (1994) shows 
that in the limit, as the number of voters grows, there are two kinds of equilibria. In both 
equilibria MV winning alternatives have the same proportion of the votes. With 
“Duvergerian” equilibria only the first loser (in M + 1-th place) receives a nonnegligible 
proportion of the votes while other losers get nothing, and with “non-Duvergerian” 
equilibria the losing alternatives all have the same positive proportion of votes. Cox (1994) 
notes that the formal findings about exact ties between winning alternatives depend on 
particular simplifying features of the model. The empirical tests he conducts of the theory 


are not quite so sharp. 


1.3. Pivot Probabilities 


A key quantity in the considered models is the probability of an individual voter casting a 
pivotal vote. Demichelis and Dhillon (2010) and Cox (1994) assume that electorate sizes 
are fixed and commonly known. This allows computation of pivotal probabilities through 
exact combinatorial computation or multinomial assumptions. Though this approach can 
yield results, computation of these probabilities are computationally demanding. In order 
to derive computationally efficient pivotal probabilities, we rely on an assumption that the 
electorate size is only known up to the distribution. In particular, we utilize the large 
Poisson games framework derived by Myerson (2000). 

Elector 7 € {1,...,N-} has incomplete knowledge about the number of voters in an 
electorate. Elector 7 receives information about the number of voters, , and assumes that 


this value is a Poisson distributed random variable with mean N: 


P(y:) ~ Pois(N) (3) 


?Cox (1994, 610) assumes that the electorate includes voters with all possible preference orders. 


Elector 7 also receives information about the proportions of voters in the electorate that 
will take certain actions: for candidate h € {1,...,H}, she learns the proportion pp, of 
voters that are going to cast a vote for h. Let ju, be the number of voters that vote for 
candidate h. Given p;, and the properties of Poisson random variables, ju; is also a random 


variable such that: 
P(tn) ~ Pois( Np) 5 Np, = Dr . (4) 


In general, we are not concerned with the raw vote counts. Rather, we want to know 
the probability distribution of vote counts relative to other candidates. Let X and Y be 


independent Poisson random variables. Then 


z 


P(X =¥ = 2) = S(e NaN) = exp(—(Ne +) (FE)? LOVIN) (5) 


where J,(-) is a Bessel function of the first kind. S(-) is the Skellam distribution (Skellam 
1946). 

Assuming Poisson vote counts, there are two situations in which 7’s vote for party h is 
pivotal: her vote pushes party h out of a tie and into a winning position or her vote pushes 
party h into a tie for a winning position. For a single winner election where the candidate 
with the most votes wins, this implies that her probability of casting a pivotal vote for 
party h (pup) is: 

Pn = > [P(Nn — Nj = ON, Nj > Nn Vm F j,h) 
j#h 
+ P(N, — N; = —1|Nn, Nj > Nm V mF j,h)] (6) 


[] 2\n — Nn > 0)P(N; — Nin > 0) 


m#zj,h 


In words, the probability that she casts a pivotal vote is the probability that she breaks or 


creates a tie times the probability that party h is at least in second place. 
Assuming that the distribution of the difference between the number of votes for two 
parties is independent of the probability that they are in the top two positions,® we can 


define pv; in terms of Skellam distributions: 


Pn = )_[S(0; Na, Nj)+S(-1; Ma, Ni) [] (os w; Ny, N, 5) [os w;N;, Nm) ) (7) 


jh m#j,h 


This value can be evaluated computationally using Skellam implementations in various 


statistical libraries. 


2 Replications 


We present examples of ABMs built on a Poisson games foundation designed to replicate 
basic results from Demichelis and Dhillon (2010), Palfrey and Rosenthal (1985) and Cox 
(1994). For Demichelis and Dhillon (2010) our goal is essentially to replicate Figure 1 in 
their paper. For Demichelis and Dhillon (2010) we stick as closely as we can to the 
assumptions they used, except for our population of electors being larger than the 
population they used to compute their Figure 1. Cox (1994) does not provide an analogous 
specific replication target, so we illustrate how ABMs can produce both Duvergerian and 
non-Durgerian equilibria. We do not restrict ourselves to utilities on the unit interval, as 
Cox (1994) assumes, but the scaling of the utilities is inconsequential for our ABM. As 
Myerson and Weber (1993) demonstrate for voting systems generally, multiple equilibria 
exist in the Cox (1994) model. We show some of the equilibria that can arise even given 


identical voter utilities. 


3In general, this statement is not true. By virtue of knowing the distribution of the vote difference for two 
parties, we also have more information about their respective rankings relative to other candidates. However, 
as N gets large and the relative differences between parties gets large, independence provides a reasonable 
approximation to the true value of puvz. Simulations and similarity of replicated results have shown that 
there are no substantive changes using this approximation over exact results. This topic will be explored 
much more in future iterations of this project. 


2.1 Strategic Abstention 


Here we focus on the simplified scenario in which the two candidates have equal support 
among electors 7 = 1,..., Ne, have the same net cost c of voting and the same initial 
probability q; = q of voting. More specifically, we initialize each run of the model 
instantiating N expected voters where half support candidate A and the other half support 
candidate B, so N4 = Ng = N/2; and we also define a value for g and cs.t. ¢q ~ U(0,1) 
and c ~ U(0,0.5). 

Let n, and ng be, respectively, the expected numbers of supporters of candidates A and 


B that would go vote if the election happened at the current iteration of the model. Then 


nj = Ds di (8) 


where i = 1,...,.N; are the electors that support candidate j. 
At each iteration t, each elector i updates her probability g; of voting, according to the 


following difference equation: 
qi(t + 1) = ait) + Ki(na(), na(t)) (9) 
To implement the function A(.) of Demichelis and Dhillon (2010) (recall (2)) we use 
Kj(na,ng) =k sign(gi(na,np) — 2c) (10) 


where constant k > 0 represents the learning speed mentioned but not specified in 
Demichelis and Dhillon (2010)—here k simply the convergence speed and run granularity 


of the model—and g;(n4, np) is the function that calculates the probabilities of elector i 


10 


being pivotal to untie or to tie the election: 


S(0,n4—1,ng)+S(—1,n4-—1,ng) if i supports A 
gi(na, NB) = (11) 
S(0,ng —1,n4) + S(—1,ng—1,n,4) if i supports B 


where S is defined in (5). S(0,...) and S(—1,...) represent, respectively, the probability of 
elector 7 being pivotal to untie the election and the probability of 7 being pivotal to tie the 
election. Given the support of the parameters of the Skellam density, we set S = 0 if na, 
na —1, np or ng — 1 is equal to zero. 

This way, as the simulation goes, the probability of any elector 7 going to vote decreases 
faster when the expected abstention is low and it decreases more slowly when the expected 
becomes higher. Up to a point in when g(.) surpasses twice the cost c of voting, in which 
case the probability of 7 voting starts to increase. The model keeps running until 
convergence is reached.* 

Figure 1 shows the final probability q of voting versus the cost of voting, subsetting 
according to different levels of initial g. As it can be clearly seen, when initial g = 0 no 
elector ever increases such probability and nobody is expected to go vote. For all other 
levels of initial q, all electors are expected to vote when c = 0, many electors are expected 
to vote when c is very small, but the expected turnout rapidly decreases as c increases, to 
the point of quickly approaching zero turnout. This result replicates the main equilibrium 
condition in Demichelis and Dhillon (2010)’s Figures 1 and 2, as well as the stable lower 


bound equilibrium. 
*** Figure 1 about here *** 


Figure 1 in Demichelis and Dhillon (2010) shows an equilibrium with full turnout for all 


4Observed convergence in this model has been of two types, which we call hard convergence and bouncing, 
hence we needed two convergence criteria to stop the simulation. A given run was considered to hard-converge 
if the values of n4 and ng remained the same for 100 iterations. Otherwise, a run could be stopped due to 
achieving stable bouncing if n4 and ng kept alternating between two exact same values every two iterations 
for 100 iterations. In both cases, since we are operating with floating point arithmetic, we consider equality 
over iterations as being equality up to &/100. 
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values 0 < ¢ < 1/2, but we do not reproduce this equilibrium. We find full 

turnout—q = 1—is a stable equilibrium for small values of c, but for larger values of c 
using an initial value of gq = 1 leads to the low turnout equilibria. Such a result is not 
surprising. As Palfrey and Rosenthal (1985, 71-73) prove, the wide range of full turnout 
equilibria does not exist for electorates as large as we are using. Whenever n, or ng is not 
small, the Skellam distribution function returns small values so that g(.) < 2+ c¢, which 
forces q to decrease. 

Of course, the fact that with all initial g greater than zero the relation between final q 
and c seems identical does not mean that the process to reach such similar equilibrium is of 
identical speed. As it can be seen in Figure 2, the distribution of in which iteration the 
model did converge follows a clear pattern. The greater the initial g, the longer the model 


takes to achieve the near zero turnout that it achieves as c gets further from zero. 


*** Figure 2 about here *** 


2.2 Wasted Vote Logic 


This paper’s replication of the Cox (1994) wasted vote model functions as follows. First we 
specify a number of candidates and a number of electors. For each elector we generate one 
utility value per candidate, u;. These utilities are generated according a user-specified 
distribution, and each elector has a utility value for every candidate. These utilities are 
assumed to be von Neumann-Morgenstern consistent and indicative of at least a weak 
preference ordering over the candidates.°® Voters then start by choosing any candidate 
other than their least preferred candidate to be their initial choice. We let voters “state” 
any nondominated initial choice to capture effects of variation in the initial beliefs that we 
have not yet modeled explicitly. Once they make this initial choice, they use public 


information about the distribution of all voters’ current choices to make a strategically 


>Given the double-precision pseudorandom numbers we use to generate utilities, no two utility values 
should be the same. 
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optimal choice, v;, according to 


uj(h) = dL Tha (uilh) — ui(Z)) 


vu; = argmax|[u;(h)] 
h 


where T},; is the probability that an elector casts a pivotal vote for h against j. Using (7), 
we set T),,; = pvp, with N, being the current “stated” vote counts for candidates 
| mel Cann 

This procedure is iterative: each elector simultaneously updates her choice to be 
strategically optimal conditional upon a) her utilities and b) every other electors’ current 
choice. But after each step the state of the world has changed, so the strategic decisions 
that were conditioned on other electors’ previous choices are no longer strategically optimal. 
So in the next iteration of the model every elector simultaneously updates her choices given 
each other electors’ current choices. The electors continue to update in this manner until 
no electors’ strategic choices change for an arbitrary user-selected number of iterations. At 
this point—and not before this point—electors’ choices are in a strategic equilibrium. 

Between any two runs of the model, there are therefore two elements of randomness. 
The central source of randomness is the distribution by which electors’ utilities are 
initialized. The second source of randomness is the arbitrary selection of voters’ initial 
candidate choice. Different results can be obtained by changing the distribution of electors’ 
utilities, or by holding those utilities fixed while rerunning the model with different initial 
choices. The first option is a completely different initial setup, which generates electors that 
have meaningfully different preference structures. The second option is a perturbation of 
the same initial electors and candidates, changing only the electors’ initial selection. A run 
of the model with the arbitrary initial selection varied should yield an equilibrium in the set 
of equilibria which are obtained by perturbing a model slightly away from its initial setup. 


With straightforward specifications of the model, our replication yields two different 
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outcomes which were derived by Cox (1994): the 2-candidate or M + 1 result (in which one 
candidate wins, another candidate gets second place, and every other candidate falls to 
near-0 votes), and a non-Duvergerian result (in which one candidate wins and all of the 
other candidates tie for last).© We demonstrate a variety of these results by rerunning the 
model with several different utility-generating distributions, and checking which equilibria 
can be obtained by perturbing the initial choice of favored candidate. 

Some equilibrium vote totals are shown in Figure 3. Each row of the figure is a different 
distribution for generating sincere utilities, and within each row are several plots which 
show the final vote totals of 4 competing candidates at the different equilibria which can be 
obtained by varying the initial selection of candidates without altering electors’ utilities. 
Generating sincere utilities according to a power law (exponent 2), logistic (uw = 2,s = 1), 
or beta distribution (a = 0.01, 6 = 0.01) typically produces an M + 1 result. Using a 
standard normal distribution or a uniform distribution typically produces the 
non-Duvergerian result, but it is also possible to obtain the M + 1 result for some 
initializations.’ It is generally possible to obtain multiple equilibria by varying electors’ 
initial candidate choices without altering utilities. These equilibria can be quite different; 
the first row shows that, by perturbing the initial choice of a group of electors that 
otherwise does not change, it is possible to make any candidate win. Several different 
equilibria can be obtained in nearly any run of the model using only initial perturbations. 
Given Myerson and Weber (1993)’s results about the existence of many equilibria in almost 


any voting system, such results are what we expect. 


*** Figure 3 about here *** 


®This code is in an Alpha stage; the results that Cox (1994) predicts make up the overwhelming results 
for simple parameter choices, but they are not fully robust to sweeps over the random parameters. The 
conditions in which the model produces unexpected results are themselves interesting, and are one of the 
top priorities for future work in this project. 

"In our replication the disfavored candidates typically fall to between 0% and 5% of the vote. This is 
presumably because this is a finite approximation of an asymptotic phenomenon, and in the very large-elector 
case those trailing candidates would approach 0%. 
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3 Discussion 


Using agents in a many-equation dynamical system and a Poisson game information 
structure, we replicate basic equilibrum results from Demichelis and Dhillon (2010), Palfrey 
and Rosenthal (1985) and Cox (1994). Symmetric mixed equilibria that Demichelis and 
Dhillon (2010) show are dynamically stable given their learning model based on a 
differential equation we find are dynamically stable using, effectively, difference equations. 
We find that both Duvergerian and non-Duvergerian equilibria such as Cox (1994) studies 
are dynamically stable in our system. Even though we do not yet explicitly model the 
beliefs that can support multiple equilibria in the wasted vote model, we find multiple 
equilibria by using different initial values for voters’ vote choice probabilities (compare 
Bouton and Gratton 2015). 

Finding these equilibria requires careful attention to numerical details. Using the 


Poisson game approach with Skellam distributions, pivotal probabilities are extremely 


small with large electorates. The electorate size we use here—N = 1000—is not large 
compared to the country-scale electorates that will ultimately concern us. Use of extended 
precision numerical libraries is crucial. In our simulations we still encounter some 
anomalous behavior that we suspect may trace to numerical imprecision. It is also possible 
that we will need to abandon our approach using difference equations—that is, each agent 
updates at discrete time steps—and move to a differential or difference-differential scenario. 
Perhaps agent update continuously, but information about expected election outcomes 
updates discretely. 

We do not yet incorporate beliefs explicitly into the models. Agents also have fixed 
preferences. In the current implementations choices are updated without reference to any 
prior beliefs—other than whatever beliefs are implicit in the specification of the current 
state—based on public information about the currently expected outcome of the election. 
We plan to augment the models to make agents explicitly Bayesian, and we will investigate 


what happens when information is more limited. For example, agents may learn only about 
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the intentions of their neighbors, in various senses of “neighbor.” Or agents may be only 
adaptively rational. 

We encountered apparent anomalous behavior in scenarios, not reported here, where we 
moved far away from the conditions specified in the formal theory papers. For instance, 
we’ve run iterations of the strategic abstention model with initial vote probabilities q; 
drawn using a uniform distribution on the unit interval. Demichelis and Dhillon (2010) say 
their results hold for mildly dispersed initial g; values. When initial q is uniformly 
distributed on (0,1), phenomena very different from their equilibria often arise. A next 
thing we will do is explore how much initial g; can vary across individuals before their 
equilibrium results break. We will also investigate variation over individuals in net costs. 

After similarly exploring the boundaries of weakening the assumptions in the wasted 
vote simulation, we expect to combine the strategic abstention and wasted vote models. 
Such a combination should produce new results going beyond what has been possible using 


paper-and-pencil formal theory tools. After that, we’ll see where it goes. 
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Figure 1: Probability of voting (q) versus cost of voting (c) 
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Figure 2: Iteration when model run converged 
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Figure 3: Wasted Vote Model with Varying Utilitiess 
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